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I Abstract 

> ■ 

OO ■ The spherically symmetric volume operator is discussed and all its eigenstates and 

eigenvalues are computed. Even though the operator is more complicated than its 
homogeneous analog, the spectra are related in the sense that the larger spherically 
symmetric volume spectrum adds fine structure to the homogeneous spectrum. The 
formulas of this paper complete the derivation of an explicit calculus for spherically 
symmetric models which is needed for future physical investigations. 
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tfo 1 Introduction 

In complicated theories such as general relativity or its possible quantization, loop quantum 
gravity PQ, much information can usually be gleaned by studying simpler situations in 
symmetric contexts. Classically one thus obtains special exact solutions, which are relevant 
if perturbations around them are stable. In the quantum theory, symmetric models always 
present an approximation, whose relevance can be investigated by comparison with less 
symmetric situations. Examples for exact symmetric solutions were found early on, with 
applications in cosmology (isotropic space, [2J Ej) and black holes (spherical symmetry, 
|3]). Quantum gravity in the Wheeler-DeWitt form has, similarly, first been studied in 
isotropic models |5j followed later by inhomogeneous ones 

In loop quantum gravity there is a systematic procedure to derive symmetric models 
and to relate them to the full theory [7J |S|. Also here, applications were first done in 
the simplest, isotropic context followed by anisotropic models. Applications include a 
general non-singular behavior in those models [3 EH EH E] , new information about initial 
conditions [13 , and several phenomenological scenarios |IH ITo] ITrj]. 

Spherical symmetry presents an inhomogeneous model with infinitely many kinemat- 
ical degrees of freedom, but also interesting physical applications for black holes. At the 
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same time, results from homogeneous models can be checked by comparing with a more 
complicated, inhomogeneous context. In many respects, spherical symmetry yields the 
simplest inhomogeneous model. Classically, there are only finitely many physical degrees 
of freedom and the system is exactly soluble |TZl HE] . The loop quantization also presents 
some features that can be expected to simplify calculations, such as the existence of a flux 
representation JH] which was helpful in homogeneous models. We will see here that also 
the volume spectrum can be obtained explicitly which is essential for calculations involving, 
e.g., the Hamiltonian constraint. However, the volume will already be more complicated 
than in homogeneous models, which will appear as special cases. 

One of the aims of studying different symmetric models is to weaken the degree of 
symmetry step by step, in order to check the robustness of results. Here, our focus will 
be on the effect of level splitting (of volume eigenvalues; for the area operator see j20j ) 
which is a consequence of weakening the symmetry and familiar from atomic or molecular 
spectroscopy. A difference in the case of quantum geometry is that the symmetry can only 
be reduced in discrete steps rather than by a continuous deformation. Mathematically 
this is implied by the fact that symmetric states are not ordinary states but distributions 
in a less symmetric model. Yet, the effect is visible, e.g. when comparing isotropy with 
anisotropy. It will be seen similarly, but in a much more finely structured way, when 
comparing homogeneity with inhomogeneity, even when only a single vertex is considered. 

This effect illustrates properties of the volume and may also shed light on the full 
case. Further applications, which will be studied elsewhere, include, e.g., the Hamilto- 
nian constraint, observables, horizons in quantum geometry. The calculations will be more 
cumbersome than in homogeneous models since the volume operator has eigenstates dif- 
ferent from those of flux operators, and the transformation between the eigenbases is quite 
involved. We will present here some of the formulas which can be helpful for explicit 
investigations. 

2 Spherical symmetry 

The spherically symmetric model from the point of view of quantum geometry has been 
described in jTU], where its states and basic operators have been derived from a symmetry 
reduction. Here we collect the information relevant for studying the volume operator. 

2.1 Classical reduction 

Spherically symmetric connections and densitized triads on a manifold X = B x S 2 can be 
written as 

A = A x (x)A 3 dr + (A 1 (x)A 1 + A 2 (x)A 2 )dfi + (A 1 (x)A 2 - A 2 (x)A 1 ) sin Mip + A 3 cos $d(p (1) 
and 

E = E x (x)A 3 sm$— + (E 1 (x)A 1 + E 2 {x)A 2 ) sintf— + {E\x)A 2 - E 2 (x)A 1 ) — (2) 
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with real functions A x , A±, A 2 , E x , E 1 and E 2 on a one-dimensional, radial manifold B 
with coordinate x (in addition to polar coordinates d and tp). The su(2)-matrices Aj are 
constant and identical to n = (with Pauli matrices 07) or a rigid rotation thereof. 

The symplectic structure is given by 

Q B = — ^- / dx(dA x A dE x + 2dA! A dE 1 + 2dA 2 A dE 2 ) (3) 
2 7^ Jb 

with the gravitational constant G and the Barbero-Immirzi parameter 7. 

The reduced theory has residual U(l) gauge transformations corresponding to rotations 
around A 3 and generated by the Gauss constraint 

G[X] = [ dx\(E x> + 2A 1 E 2 - 2A 2 E X ) « . (4) 
Jb 

The components (Ai, A 2 ) and (E l ,E 2 ) transform in the defining SO(2)-representation, 
such that the components 



Mx) := vW + W- (5) 
:= v^W+W (6) 

are gauge invariant. The corresponding internal directions 

A£(x) := (Ai(x)A 2 - A 2 (x)Ai)/A v (x) , (7) 
Aj(x) := (E 1 (x)A 2 -E 2 (x)A l )/E^(x) (8) 

are rotated simultaneously such that 

cosa(a;) := A£(x) • Aj(x) (9) 

is also gauge invariant. Parameterizing 

A£(x) =: A lC os/3(x) +A 2 sin/3(x), (10) 
Aj(x) =: Aicos(a(a;) + + A 2 sin(a(:r) + (3(x)) (11) 

introduces the angle (5{x) which is pure gauge. 

In these variables which take into account the gauge structure we have the symplectic 
form 

Q B = I dx (dA x A dE x + 2dA v A d(P^ cos a) + 2dj3 A d(A ip E Lp sin a)) 

2 7^ Jb 

= ^— [ dx (dA x A dE x + dA v A dP 1 ^ + d(3 A dP' 3 ) (12) 
Jb 

with momenta 

P^x) := 2P^(a;) cosa(:r) (13) 
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conjugate to A v and 

pP(x) := 2A v {x)E ip {x) sin a(x) = A^(x)P^(x) tan a(x) (14) 
conjugate to (3. The Gauss constraint then takes the form 

G[X] = [ dx\(E xf -P^xiO (15) 

J B 

which is easily solved by P 13 = E xl while the function A x + f3' is manifestly gauge invariant. 

Thus, in this model flux variables (i.e., momenta of the connection), are not identical 
to triad components. In particular, some metric components in 

ds 2 = E x {x)~ 1 E^(x) 2 dx 2 + E x (x)(d$ 2 + sin^dy? 2 ) , (16) 

and in particular the volume V(I) = 4ir j x dx^E^E^ of the shell IxS 2 cE given by an 
interval X C B, are rather complicated functions of the momenta as well as of connection 
components. 

2.2 Loop quantization 

States in the connection representation depend on the components A x , A v and (3. The loop 
representation is based on cylindrical states with respect to graphs g in the one-dimensional 
radial manifold B, consisting of a collection of non-overlapping edges e whose endpoints 
define the vertex set V(g). We assume a given orientation of B and choose the same one 
for all edges. Each vertex v G V(g) then has one outgoing edge e + (v ) and one incoming 
edge e~{v ). An orthonormal basis of states is given by spin network states 




in the gauge invariant form, with labels k e G Z for U( ^-representations on all edges e C B 
and fi v G M for all vertices v G B. Thus, in addition to ordinary holonomies of A x along 
edges there are vertices with point holonomies exp(i//„A ¥ ,(„0 which can be thought of as 
representing orbital edges in the symmetry orbits S 2 . Gauge invariance requires only that 
the differences k e +( v \ — k e -t v ) for any two edges e + {v) and e~(y) meeting in a common 
vertex v must be even. Otherwise the edge and vertex labels are free. 

Connection components are quantized via holonomies and act directly by multiplication. 
For instance, looking only at a single vertex v and suppressing all labels but fi v , we obtain 

sin(M^ = -§*(T^-T^) (18) 
008(^)7^ = KW + ^-m) (19) 
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as basic operators representing the transversal connection component which takes values 
in the Bohr compactification Meohr of the real line. Flux components act as 

E*(x)T gA , = ^ W ^- W r^ (20) 

pH>T 9,k,iX = ^2 ^ vT 9,k^ ( 21 ) 

veinv(g) 

P^Tgfi^ = ^ k v T g ,k,tM (22) 

v<=XnV(g) 

with fc„ := — fe e -( w )) G Z, which shows that they all have discrete spectra (normal- 

izable eigenstates). 

The momenta P^ and P^ transform as scalars of density weight one and thus have to be 
integrated over intervals X C B in order to yield well-defined flux operators. Alternatively, 
they can be regarded as operator valued distributions 

P' p (x)T gAll = -^- ^ l^vS(x,v)T gAlM . 

71 v<EV(g) 



3 Volume operator 

Classically, the volume of a spherically symmetric region between two spherical shells 
located at points X\, x% 6 B is given by 

V(T) = 4vr J Ax^E^E* (23) 

where X = [xi,x 2 ]. Since \E X \ is a non-negative self-adjoint operator, we can directly use 
it in this expression. The triad component E v , on the other hand, is more complicated 
to quantize since it is related to the basic variables by (|T3~|) which also involves connection 
components through a. Using (JT3j) and (|14j). we obtain 

= -J— = ±P*Vl + tan 2 a = U(Pv) 2 + A- 2 (P?) 2 . (24) 
2 cos a V r 

In order to quantize E v and diagonalize the resulting operator it suffices to look at the 
action on eigenstates of E x and thus P 13 since this does not restrict the /^-dependence 
of states. In the expression for E^ we can then regard P@ as a constant / y£pk/4Tr and 
use basic operators for P^ and A^. Thus, (E^) 2 has the form of a classical mechanics 
Hamiltonian with potential k 2 j A 2 ^. Taking this Hamiltonian at face value suggests that 
the volume spectrum is continuous since the potential does not have an absolute minimum 
at finite values. In fact, since E v depends on fluxes as well as connection components 
there is no a priori reason to expect a discrete spectrum. However, a continuous spectrum 
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would be difficult to reconcile with the full volume operator (2JJ as well as with results in 
homogeneous models ITU} ITTj. 

This problem is resolved easily after noticing that there is no operator for directly 
in a loop representation, or on the Bohr Hilbert space in this case (see the discussion of 
in the isotropic context). Then neither the quantization of A~ 2 as a multiplication 
operator on a Schrodinger Hilbert space, as understood in the above argument, is available. 
In fact, the discreteness of geometric spectra depends crucially on characteristic features of 
the loop representation. Following a loop quantization of and checking if the resulting 
spectrum will turn out to be discrete can thus be seen as a further test of the generality 
of loop results. Since the spherically symmetric volume is a rather complicated function of 
the basic variables, but not basic itself, one cannot expect a unique (up to smaller choices) 
quantization. But the sensitivity of general properties like the discreteness of its spectrum 
and its behavior can easily be studied with different versions and compared with the full 
theory as well as other models. 

In a loop quantization only holonomies of the connection are represented as well-defined 
operators on the Hilbert space, not connection components directly. This is also true for 
connection components acting on a Bohr Hilbert space as in (jl8j) . (fT§|) . Since the states (fTTj) 
are orthonormal, the /z-derivative at fi = of the operators sin(/iy4 v ,) and cos(/xA ¥ ,) does not 
exist such that an operator for cannot be derived in this way. Instead, we have to use 
(118)1 and (j!9)) with non-zero fi as our basic operators in every expression where A v appears 
and needs to be quantized. This in general is true for the Hamiltonian constraint which 
contains curvature components of the connection. A quantization can then be constructed 
by using holonomies of the connection and in a homogeneous model or for orbital 
components of an invariant connection point holonomies j2Hj- Since in the spherically 
symmetric model the volume is more complicated in terms of basic variables, containing 
in particular connection components, we have to use the same technique here. We will 
exploit this fact in order to study consequences of replacing connection components by 
holonomies by means of comparing the spherically symmetric volume spectrum with that 
of homogeneous models where such a replacement is not necessary for the volume operator. 

The connection components contain information about intrinsic curvature (via the spin 
connection) and extrinsic curvature. In a semiclassical regime, which we need to take 
into account to ensure that the resulting volume operator will have the correct classical 
limit, extrinsic curvature is small but the spin connection may not be small when we 
restrict ourselves to spherical coordinates which just leaves the freedom to transform x 
(under arbitrary coordinate transformations the spin connection can locally always be 
made as small as desired) . For our purposes it is nevertheless sufficient to assume that the 
connection components, or at least A v , are small semiclassically since the spectrum will 
be invariant under shifting A v by a constant in the operator used later. 

We thus replace A v by sin(5^4 ¥ ,)/5 for some 5 G M\{0}. In the limit 5 — > we would 
obtain the classical expression exactly, and for 5 ^ we have small corrections in regions 
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where A v is small. This leads us to consider the Hamiltonian operator 




2 



H k>s = (47rP7 7 4) 2 + 5 2 fc 2 sin(<L4„) 



-2 



(25) 



for integer k and non-zero real 5. It will later be used to quantize E v and the volume. At 
this point we can already see that replacing by the sine, which we are forced to do in 
the loop representation, leads to a discrete spectrum since the potential is now unbounded 
above at both sides. This fact does not depend on the explicit realization of the replacement 
and is thus robust to quantization ambiguities. In order to test the replacement we will 
later use more detailed properties of the spectrum and compare it with volume spectra in 
other models. 

Before we continue we have to discuss the Hilbert space on which the operator ()25|) 
acts. Since it suffices to consider a single vertex and we already used the eigenvalue k of P 13 
it acts on a single copy of the Bohr Hilbert space with A^ G ffiBohr- Thus, the Hamiltonian 
is not yet of the usual quantum mechanical form since so far it is not represented on 
the Schrodinger Hilbert space. However, it is clear from (fTK|) that H^^ leaves invariant 
any subspace spanned by states exp(i(n<5 + e)A Lp ) 1 n G Z with fixed e G [0, 5). Each such 
subspace is isomorphic to the Schrodinger Hilbert space on a circle where for e^Owe have 
to use a non-trivial line bundle for the states. That the inner product on the subspace 
spanned by states exp(i(n5 + ejA^), n G Z is identical to that for states on a circle follows 
from 



using the Bohr measure and substituting y = 5A !P /2tt and m = 5T/2tt. Now we have 
an ordinary quantum mechanical system with Hamiltonian (|23j) whose spectrum can be 
determined with standard methods (from this perspective, e plays the role of a #-angle). 

3.1 Spectrum of H 

Using P v oc d/dAp and substituting 6A V for A v , the operator family (|23|) labeled by 5 
satisfies the identity Hk,8 = S 2 Hk,i such that it is enough to determine the spectrum for 
5=1. Moreover, since the sequence of states running through P^-eigenvalue zero is most 
important, we will focus on the case e = 0, for which states are periodic functions on 
a circle. There are then standard techniques (see, e.g., [251, P- ^32) to determine the 
spectrum of the operator 




T^oo 2T J 



lim I e - t(ni - n2)SA *dA, 




5„ 



ni,n 2 




d£ 2 



+ 
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where for later convenience we introduce p := <J k 2 + | > which for our case of integer k 
is always non-integer. 

If we know an eigenstate ip? with eigenvalue A„ labeled by some parameter n, i.e. 
Hpipn = ^nV'ni we can fi n d eigenstates of H p ±i by applying the operators 

a±:=±A + (p_!) co t£ (26) 



to the state ip?. (For our values of p, p ± 1 will not be of the form yk 2 + ^ with integer k. 

Nevertheless, these values can be used in intermediate steps to construct the eigenstates.) 
This follows from the identities 

a-a+ = H p -(p-\y (27) 
K+rK+i = H p -(p + \f (28) 

which imply 

H P -ia + p = {a+a p + (p - |) 2 )a+ = a+ (a~ a+ + (p - |) 2 ) = a+H p 

and similarly 

Thus, if xjjp is an eigenstate of H p with eigenvalue A n , then a p ipp is an eigenstate of li p -\ 

and aL.iV'n an eigenstate of H p +\ with the same eigenvalue A n . 
In order to normalize the new eigenstates we compute 

j \% + lV n \ 2 ^ = f V^fijhl WW = ( A « - ^ + / ^n\ 2 ^ ( 29 ) 

and 

'|a^| 2 de = J V:a' p a + p Vn^ = {K-{p~\f) J \V n M ■ (30) 

Since all eigenvalues of H p are larger than the potential minimum p 2 — \ = (p+|)(p— |), we 
must have started with an eigenvalue A n > (p — |) 2 . When applying a~ several times, we 
get eigenstates to H p with higher and higher values for p, and since A n stays the same in this 
process and the integrals in (|2*5|l are non-negative, (|2*5|l is consistent only if A n = (p+n + |) 2 
and n G No a non-negative integer. The state with the highest allowed p at fixed A n , i.e. 
for n = 0, is then annihilated by a~ +l as a consequence of 



which has normalized solution 



'fr(2p + 2) 

" 2^r(p + i) s v ; 



Normalized eigenstates with values n > are obtained by applying a+ and using (|30|) 
with X n - - i) 2 = (2p + n)(n + 1) such that = [(2p + n)(n + l)]" 1/2 a+^. The 
eigenstate ^ can thus be obtained from the known ^>Q +n as 



r(2p + n + l) +r , 
' " V 1(2/, • 2n ■ 1) Vl 

Using 



iL + (p- I) cote = (sine) 1/2 -^(sine)"- 1/2 = -(ffln0 s ^'- r ^(idn0'^ 
d£ d£ dcos£ 



such that 



this yields all eigenstates as 



(p + n+ i)r(2p + n+ 1) H n 

= (-If -v ^siii^eTT^Trsin^e (32) 

2"+™r(p + n + l)^! S (dcosO" 



of .f/p with eigenvalue \ n = (p + n + 1/2) 2 , n e N . If p = 1/2, i.e. k = 0, the potential 
vanishes and the poles at x = 0, it disappear. In this case boundary conditions leave more 
freedom since states do not need to vanish at x = 0, tt which results in the additional 
eigenvalue A = for the case p = 1/2. 

These eigenfunctions can be expressed in terms of Legendre functions 1 P£(z) (see Ap- 
pendix) 

Vni.0 = \J{P + n + \)Y{2p + n + 1)/ n\ y/\ sin£| P^cosf) (33) 

such that the sum of degree and order, v + p = n, is a non-negative integer and the order 
p = — p is negative. 



3.2 Volume spectrum 

We now have to construct the volume operator from the basic building blocks E x and 
E v = -^-yH. In particular, expressions like (P^) 2 and (P' 3 ) 2 in H have to be regularized 
upon quantization since P v and P 13 are operator valued distributions. 

We do that by decomposing a given interval X, whose volume V = Air J x dxy/\E x \E ip 
we want to quantize, into sub intervals Zj = [xi — e/2, X{ + e/2] each of size e: X =\} i X i . 

1 It is amusing to note that such as in rotationally symmetric systems in quantum mechanics also in 
spherically symmetric quantum geometry Legendre functions play an important role (with the difference 
that here the labels are always non- integer). 



9 



The volume can then be written as 



V(J) = 4tt / dxy/\E?\E' f ' = 4vr V / dar^/fE*^ 

i i 

where now only squares of fluxes occur which can immediately be quantized. We now 
choose the subdivision of X fine enough such that each Tj contains at most one vertex of 
the spin network state we act on, and immediately obtain the volume operator 

V£) = ^ E \/\Hv)\H{v) (34) 

V 

where only vertices of graphs appear since otherwise P^ and P 13 , and thus H(v) = 
—d 2 /dA ip (v) 2 - (sin A^v))" 2 d 2 /d(3(v) 2 , would annihilate the state. 
Using the spectra of E x and H, we find the volume spectrum as 

V k , n = (| 7 4) 3/2 (47T)" 1/2 J2 \l\\ k eHv)+ke~{ V )\ (n v + \ (l + ^/l + ( k e+(v) - fc e - („) ) 2 ) ) 

v 

(35) 

where fc e +(«) is the label of the edge outgoing from v, k e -^ the label of the edge incoming 
to v, and n v a label of the vertex with n v G N if & e +(«) 7^ &e-(u) an d n„ e N U { — 1} 
if k P+ (v) — k e -( v y Moreover, for the spectrum on gauge invariant states the difference 
k e +{v) — k e -M must be even [TH], and thus also the sum k e +f v ) + k e -r v )- 

The A^-dependence of volume eigenstates is now via an if-eigenstate xp^JA^), where 

p(v) = |v/l + (k e +(v) — ^e-(i))) 2 using the eigenvalue §(& e +( v ) — fc e -(„)) of P 13 , explicitly 
given by 

S gAn = l[exp(iikJ(A x +fl')dx\ J] 4f +(K+{v) ~ K ~ {v)) \A v{v) ). (36) 
3.3 Level splitting 

As in any loop quantization, the spherically symmetric volume operator has a discrete 
spectrum. Due to the symmetry, however, the set of eigenvalues (considering also non- 
zero values for e) is the whole positive real line (see for a discussion of this issue in 
the isotropic context; discussing observables may reduce the spectrum to a discrete subset 
|26j). It is nevertheless instructive to compare the main series of the volume spectra, 
which are obtained by restricting to a separable subspace corresponding to equidistant flux 
eigenvalues containing zero, for different symmetric models. Similarly to the spectroscopy 
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1 



(i) 



(ii) 



(iii) 



Figure 1: Comparison between the main series of the isotropic volume spectrum n 3 ^ 2 , 
nGN (i), the homogeneous LRS (single rotational axis) volume spectrum m-y/n, m,n G N 
(ii) and one vertex contribution of the spherically symmetric volume spectrum (iii). 
Constant factors resulting from the symplectic structure have been ignored. 
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of atoms, one then observes that levels split when the symmetry is weakened [20], as 
illustrated in Fig. [TJ 

Comparing in particular the homogeneous LRS spectrum with a single spherically sym- 
metric vertex contribution shows that the spectra differ only in details such that the homo- 
geneous levels split only slightly. Moreover, all the levels of the LRS spectrum are preserved 
and reappear in the spherically symmetric spectrum for k e +M = k e -i v y Since the classical 
geometry in a single point of a spherically symmetric space is given by a homogeneous 
LRS metric we would in fact expect a close relationship between the quantum geometries 
in a given vertex. A spherically symmetric vertex is closest to a homogeneous state if 
k e +(v) = fce-(«) ( a homogeneous state (23 HI] would be built with point holonomies also of 
A x , which can be seen as holonomies for closed edges leaving and entering in the vertex 
which requires equal spin numbers) in which case the spectra are identical. If k e +( v ) 7^ & e -(«); 
inhomogeneity becomes noticeable and the levels split. One can view k e +i v \ — k e -( v \ G 2Z 
as a new quantum number of spherically symmetric vertices characterizing the degree of 
inhomogeneity. The larger k e +^ — k e -r v \, the more strongly the levels split. Even though 
there is no smooth transition between a homogeneous LRS state and a spherically sym- 
metric one since a more symmetric state is distributional in the less symmetric model, one 
can visualize the splitting by introducing a continuous parameter h G [0,1]. With this 
parameter, we have a family of spectra 

<Y,n = V ^ + + M (" + \ (l + y/l + h\k + -k.y)) . 

For h = we get the homogeneous spectrum, while h = 1 gives all levels of the spherically 
symmetric vertex contribution. Values in between do not correspond to physical models but 
are introduced for the sake of illustration. This clearly shows that levels with k + — k_ = 
do not split, whereas all other levels split approximately linearly and the stronger the larger 
k + — is. Fig. |21 shows the splitting of all levels which will be smaller than five after the 
splitting is completed, from which one can recognize the lower part of Fig. ^ 

Thus, even though the spherically symmetric volume operator at first sight looks very 
different from the homogeneous one (where fluxes and triad components are identical and 
holonomies are not used for the operator), its single vertex spectrum displays only slight 
traces of the inhomogeneity. As discussed, there are several ambiguities affecting the 
spherically symmetric operator, mainly coming from the replacement of A v . If we choose a 
value for 8 different from 5 = 1 as used above, the spectrum of H is rescaled by S 2 and that 
of V by 5. For this value, we would be considering a different series of states given by e mSA,fi , 
which would similarly define a different series of states in the LRS Hilbert space. On those 
states, also the LRS volume operator is rescaled by exactly the same value 5 such that 
the relationship between the spectra is preserved. If we choose a different function other 
than sm5A v , on the other hand, the spherically symmetric spectrum would change 
without a corresponding change in the LRS spectrum. This would have an effect on the 
exact position of the new levels in the spherically symmetric spectrum compared to the 
LRS one. 

Comparing all three cases of Fig. ^ shows that the difference between the isotropic 
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spectrum and the homogeneous LRS spectrum is much more pronounced than that between 
the homogeneous LRS spectrum and the single vertex spherically symmetric spectrum. 
When breaking isotropy, many more lines emerge which cannot even be attributed to 
individual level splitting, and a spectrum with increasing level distance turns into one with 
decreasing distance. When breaking the homogeneity of the LRS model, on the other hand, 
the single vertex spectrum results from small level splitting, which can be considered as 
just fine structure of the overall spectrum. This may be taken as an indication that the 
transition from homogeneous to inhomogeneous models, though technically very involved, 
is less severe in geometrical and physical consequences than the transition from isotropy to 
anisotropy. Since physical results such as the non-singular behavior [Oj have been extended 
from the isotropic context to the anisotropic one [TTJ ^] , an extension to inhomogeneous 
models may be feasible, too. 

Note, however, that we used only a single vertex contribution of the spherically sym- 
metric volume spectrum for this comparison. Since the spectrum (J35j) is a vertex sum over 
non-integer contributions, it will become dense much more rapidly at large eigenvalues 
than just a single vertex contribution. 

3.4 Triad spectrum 

Eq. shows that the eigenvalues of the triad component E v in a single vertex are given 
by the square root of those of the operator H, 

v 

Unlike for the triad component E x or flux operators, this spectrum is not equidistant since 
vertex contributions are in general irrational. Thus, for large eigenvalues the behavior is 
very different and the spectrum crowds rapidly if several vertices are considered. This 
observation may be of relevance for the semiclassical limit. 

It is possible to choose a quantization resulting in an equidistant spectrum even for the 
triad component E v . To do so, we replace (P 13 ) 2 in by (P 13 ) 2 — 7 2 £p/647r 2 such that 
the label p used in the derivation of the iPspectrum would now be identical to the integer 
k labeling P^-eigenstates. Such a replacement would not spoil the classical limit since the 
added constant vanishes when H — > 0. Vertex sums such as in the P^-spectrum would then 
be over half-integers, which would not crowd even for many vertices. Moreover, a vertex 
contribution to the spherically symmetric volume spectrum would be identical to the LRS 
spectrum such that all fine structure would disappear. 

However, this replacement only works on states with non-zero P^-eigenvalue, for oth- 
erwise the potential would be negative and unbounded from below when k = 0. One can 
avoid this by choosing (P' 3 ) 2 — 7 2 £p/647r 2 on P^-eigenstates with non-zero eigenvalue, and 
the previous prescription (P 13 ) 2 on the zero eigenstate. Since this special treatment of the 
zero eigenstate seems artificial, we do not use the resulting operator here and regard the 
non-equidistant triad spectrum and the related fine structure of the volume spectrum as a 
robust consequence of the loop quantization. 
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4 Operators in the volume eigenbasis 



We now have the full eigenbasis of the volume operator at our disposal which can be used 
for further calculations. As another consequence we have seen that the triad operators E x 
and Jj E v commute with each other such that a triad eigenbasis exists (given by the volume 
eigenstates). This may be useful for investigations of the Hamiltonian constraint just as in 
the homogeneous case fOJEJE]- The difference to the homogeneous case is that the triad 
representation is not identical to the flux representation which was derived for spherical 
symmetry in j!9j . For particular calculations one then has to see which representation is 
easier to use, and maybe has to transform between the triad and the flux representation. 
We will use the explicit form of the triad eigenstates together with identities for Legendre 
functions (see, e.g., 28] ) in order to write down this transformation and to express some 
operators in both the flux and the triad representation. 



4.1 Volume eigenstates in the flux representation 



For the Legendre functions we have the trigonometric expansion (|A.9|) which gives the 
A^-dependence of volume eigenstates for fixed k + , fc_ and n, and thus fixed p and n as 

r n {A v ) = 2 1 ^7r- 1 / 2 y / n!(p + n + \)T(2p + n + l)(sin^)^ +1 / 2 

The expansion in flux eigenstates is then obtained by using 

isin/y-^/ 2 = £(-!)* S k lP ) ™s 2k A, (37) 



k 

k=0 



and 



as well as 



cos jAp sin kA v = i(sin(j + fc)Ap — sin(j — /c)^ 



after which collecting terms with equal frequency leads to the expansion of ip^A^) in terms 
of flux eigenstates. 



4.2 Holonomies and fluxes in the triad representation 

Holonomies cosA^ and sinA^ and the flux d/dA^ have simple actions (fTT)|) . (fTHjl and (|2*T]) 
in the flux representation, but are more complicated when acting on volume eigenstates. 
Still, in some cases one can make use of some identities satisfied by Legendre functions 
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(which for convenience are collected in the Appendix) in order to get explicit expressions 
for the action of holonomy and flux operators in the volume eigenbasis. 

The simplest operator will be cos A v since we can use (IA.1J) and obtain a simple recur- 
sion where p is fixed: 



: n(2p + n) l0 , . 

If p — |, for instance, both coefficients reduce to 1/2 and we reobtain ((19)1 since in this 
case a volume eigenstate is a flux eigenstate. 

Already sinA v is more complicated since simple recursion relations like (|A.2|) would 
change the label p by an integer, which would not be a volume eigenstate (no difference 
between different allowed values of p = a/ k 2 + 1/4 with integer k is an integer). One can 
use these functions in intermediate steps of the calculation, but in the end all have to be 
expressed in terms of the allowed volume eigenstates. Even when the sine has been used 
in intermediate steps, there may be cancellations (either directly or after using identities 
like (jA.5|) for Legendre functions) in some cases such that in the end the result will be 
expressed as a linear combination of finitely many volume eigenstates. Instead of using a 
recursion relation, we can use the expansion (J37|) and use the action of cos A v . This results 
in an infinite linear combination of volume eigenstates. 

The square sin 2 ^, on the other hand, does have a simple action since it can simply 
be written as 1 — cos 2 A v , for which we can use the action derived before. We then obtain 

-,2 



t (p + n)(p + n + 1) - 1 + p 2 
(2p + 2n + 3)(2p + 2n 



(sinM^(^) = 2 % ' ' - " VSAp) 



'n(n-l)(2p + n)(2p + n-l) 



2p + 2n-lV (2p + 2n + l)(2p + 2n-3) 



1 (n + l)(n + 2)(2p + n + l)(2p + n + 2) 



2p + 2n + 3V (2p + 2n + l)(2p+2n + 5) 



which again for p = | reduces to the flux result since all coefficients collapse to \. 

Flux operators, i.e. derivatives acting on the Legendre function, would in general also 
change p by an integer. When we multiply holonomies with a flux operator, however, there 
may be simpler expressions in some cases. For instance, (|A.8|) implies 



d^(A^) ^2p + 2n + l ( /( n +l)(2p + n + l) 



dA v 
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such that (sin A^P^ would have a simple action on volume eigenstates, even though the 
action of sin A v or P v separately is more complicated. 

Those operators are necessary to consider when composite operators such as the Hamil- 
tonian constraint are computed. One common example is that of commutators between 
holonomies and the volume operator, such that powers of sine and cosine act with the 
volume operator in between. 



4.3 Euclidean observables 

In the spherically symmetric system has been reduced explicitly in complex Ashtekar 
variables and in particular a complete set of observables has been found. We can use the 
same expressions, translated to real variables, for the Euclidean model, which are then 
given by the ADM mass 

m = {A\ + Al- l)^W\ = {Al - (39) 

evaluated at the boundary at infinity, and its conjugate 

JB y/±i -|- Si-2)v 1-^ I JB Z/ W I I 

In the original variables Ai, E 1 the expressions look rather complicated to quantize, but 
using the variables A^ and P v adapted to the loop quantization, we can quantize the 
observables along the lines used before, after replacing A v suitably by holonomies. 
We use sin A v instead of A v in ()39j) and obtain the operator 



m 



cos 2 ^(oo)vl£ x (oo)| (41) 



without ordering ambiguities. The argument oo here refers to the boundary point of a state 
corresponding to spatial infinity. For ()4()j) we choose the replacement sin A^ cos A v instead 
of a simple sine in order to preserve the conjugacy of the new expressions. Moreover, here 
we have to be careful with the ordering between A^ and P v . If we again subdivide the 
integration over B into integrations over intervals Xj of size e and with mid-points Xj, we 
obtain 

T= {2 1 K)- 1 \imY^E x {x i Y 1/2 A v {x i y 1 I dxP*(x) . 

After quantizing we will be able to remove the regulator only if there are only finitely many 
non-zero contributions to the sum. This is the case if we use the ordering as indicated above, 
with the flux to the right. For this ordering we obtain the operator 

f = Y,(\E^^ 1/2 )(smA^v) cos A^y'P^ (42) 

V 

which is not symmetric in the kinematical inner product (P^ refers to a single vertex 
contribution to (|21|) ). If we reorder before we remove the regulator, the resulting operator 
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will not be well-defined since P*'(sin A^cos A^)' 1 would have non-zero action even where 
the state has no vertex. One can, however, reorder after the regulator has been removed, 
such that we simply symmetrize the operator T. 

There are also several inverse expressions in T, which could make the operator ill- 
defined. We indicated in the notation above that we use a well-defined, even finite operator 
for l^'l -1 / 2 along the lines of The inverses of smA^ and cosA^ are obtained as mul- 
tiplication operators, which are unbounded but densely defined. Computing the operator 
T explicitly is nevertheless complicated since the flux eigenstates do not lie in the domain 
of definition of the inverse sine and cosine. 

On volume eigenstates the mass rh has a simple action since we only need to multiply 
with the cosine as in (JHHj) . For T, on the other hand, we have the flux operator which 
has a complicated expression in the volume eigenbasis, and also inverse powers of sine and 
cosine. The inverse sine can be directly applied to a volume eigenstate using the relation 
(IA.7j) . where the right hand side has changed orders of the Legendre functions such that it 
effectively corresponds to an infinite superposition of volume eigenstates, but the inverse 
cosine again does not have the eigenstates in its domain of definition. 

Looking at the algebra of the quantum observables rh and T with a quantized Hamil- 
tonian constraint and the issue of self-adjointness with respect to possible physical inner 
products can provide additional tests for the quantization scheme. In particular the fact 
that operators like the observables and also the Hamiltonian constraint have to be rewrit- 
ten in terms of holonomies makes it more non-trivial for the quantum algebra to reflect 
faithfully the classical algebra. This is in addition to factor ordering issues that play a role 
in any quantization with more complicated constraints. 

5 Discussion 

We have diagonalized the volume operator of spherically symmetric quantum geometry 
explicitly, which shows that spherical symmetry can supply a substantial amount of sim- 
plification compared to the full theory. This makes it possible to study also composite 
operators explicitly, where usually the volume operator plays a central role. Having an 
explicit diagonalization of the volume operator is one of the reasons why homogeneous 
models can be dealt with directly with by now many cosmological applications. However, 
while explicitly diagonalized, the spherically symmetric volume operator is much more 
complicated than the homogeneous or isotropic ones. Its eigenstates are not simply flux 
eigenstates but have a more complicated dependence on the connection components via 
Legendre functions. Since holonomies in general have a complicated action on volume 
eigenstates, composite operators like the Hamiltonian constraint, where commutators of 
the form h\h~ x ,V\ between holonomies and the volume operator appear, become more 
difficult to handle. Unless subtle cancellations happen in the construction of those op- 
erators, for instance making use of relations like (jA.5|) between Legendre functions, the 
Hamiltonian constraint in the triad (or flux) representation will be a difference operator 
of infinite order in the vertex labels n v (while the order in edge labels k e will be finite as 
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in homogeneous models [HOI EH E!)- The physical significance of this fact depends on the 
interpretation of solutions to the constraint and the issue of time. Compared to the full 
theory the possibility of explicit, though rather cumbersome, expressions will substantially 
improve the outlook to understand the quantum system and in particular field theoretical 
issues (the constraint algebra, semiclassical states) that appear here. At least numerical 
computations should be possible to implement in a straightforward way, while even this is 
complicated in the full theory |3*T] . 

Since the spherically symmetric volume operator is, unlike the full one as well as ho- 
mogeneous expressions, a function of the basic expressions involving also connection com- 
ponents, there are more quantization ambiguities. There are, however, no factor ordering 
ambiguities and so the consequences are rather minor compared to, e.g., quantizations for 
inverse powers of the volume in homogeneous models where the ordering is relevant, 

too. In particular the good agreement with the LRS homogeneous spectrum is reassuring, 
where quantization ambiguities just affect the fine structure of the spectrum. 
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Appendix 

A Properties of Legendre functions 

Recursion relations [2%] : 

zP H z ) = V ~^ +1 P»( Z ) + ^-P^ x {z) (A.l) 

yf#=lPS(z) = ^^(Pgfa) - P£liz)) (A.2) 
P^{z) = (z 2 - I)" 1 /* ((„ - rfzPftz) - (y + n)P^ x {z)) (A.3) 

(z 2 - 1)^^ = vzPftz) - {v + tiP^iz) (A.4) 
Combining (|A.1|) with (jA.3|) and then with (jA.2|) shows that 

" = {v-M»-vl+ l)P? +1 (z) + fi)P^(z) . (A.5) 
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From (|A.3|) together with (|A.lj) we obtain 

V^IP^Z) = ((„ - Ji + 1){V - H + 2)P^{z) -{y + fiiy + p- l)P^(z)) 

(A.6) 

which is independent of the relation (jA.2|) . We can thus combine both relations, shifting 
the label p, by two in (|A.6|) . in order to eliminate P„-x(z) on the right hand side. After 
dividing by \/ z 1 — 1, we have 

' --P£(z) = ±- {{u + + p, - \)P£i(z) - PEltz)) . (A.7) 



Most of these identities involve Legendre functions with different orders p,, in partic- 
ular those involving derivatives. One combination with fixed order can be obtained by 
computing 

d£ 2vsm£ dcos£ 



7 P(z) - (1 - z 2 ?/ A 



dP(z) 



2(1-^2)1/4 V J V / dz 

where z = cos£. Using (|A.4|) and (jA.l|) we obtain 

d(1 ~ * 2 ^ /4iy( * } = (i-zr^^+Dzp^-iu+^p^z)) 

= 1(1- z 2 )- 1/4 (iy - M + l)P^!(z) - iy + n)PS- X (z)) . (A.8) 
Trigonometric expansion 



(A.9) 

where (i/), := T(u + l)/T(u). 

B Legendre functions for negative order and integer 
sum of order and degree 

The formula 

/ i \n+m Jn+m 

1 J , (1 - z 2 ) m ' 2 ——{l - z 2 ) n = P™(z) 

2, n n\ ^gn+m v ' n \ / 

can be found in the standard literature (e.g., [2HEH]) for positive integer labels n and m. 
This would not be applicable in our case since p, = —p and v = p + n are non-integer. 
However, for all the volume eigenstates p, = —p is negative and p, + v = n is a non-negative 
integer. In this case, the same formula applies. 
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Lemma 1 Let (i < 0, and (i + v be a non-negative integer. Then 

P^(z) = - (1 - z 2 Y' 2 3~— (1 - z 2 ) v (B.10) 

Proof: We will prove (jB.10|) by induction over n — v + (i. Since [2E| 



2^r( P + 1) 



(B.11) 



the identity is true for n = 0. Let us now assume that it is true for a given n. Multiplying 
both sides of (jB.lOj) by (1 — z 2 )~^ 2 and differentiating with respect to z we obtain 

= M z(l - Z^-^P^Z) + (1 - ^)"M/2^(f) 

= itz(l - ^yl^P^z) 

+ (1 - , 2 )^ 2 - 1 (-(1/ - M + 1)PA-!^) + + i) W*)) 

= -(1 -^)- ( ^ +1)/2 P^ +1 (^) (B.12) 

where in the last two steps use was made of 

dgf(z) 1 
dz 1 - 2 2 

and 

(i/ - (i + l)P? +1 (z) - (y + n + l)^(z) = vT^i** 1 ^) , (B.14) 

respectively. Here, (jB.13|) follows directly from (jA.4|) together with (|A.1|) whereas (jB.14j) 
is identical to (|A.3jl . Multiplying equation (|B.12J) by —(1 — z 2 )^ +1 " 2 we infer that our 
identity (jB.lOJ) is also true for n+1 thus completing the proof. □ 

Again for positive integer m and n the function P™(z) (a polynomial in this case) is 
normalizable with 



((i/ + l)zP?{z) -(v-(i + (B.13) 



i 

P™(z) 2 dz 

1 



1 (n + m)! 
" '' n+ | (// - ///)! ' 



Also this formula holds true in our case: 

Lemma 2 Let (i < 0, and (i + v be a non-negative integer. Then 



1 (" + /x)l 
i t/ + i Y{y — (i + l 



P?(zfdz = .Z™ , • (B.15) 
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Proof: We use the representation (jB.lOj) and integrate by parts where boundary terms 
vanish: 



2 2l T(z/ + l) 2 J v 7 Vd^ +At 

T(i-- 2 )")((i-- 2 )"CT(i-^r 



2 2i T(z/ + l) 2 J Vd^^- 1 /V d^+^+ 

^(i-- 2 ) ( - 1)/2 (^=i(i-- 2 r) 

dz 2 V Vd^ +M 

-2z(l - zH^yz^-n - Z 2 Y 

+{H - 1)(1 - (ai - 3)^/(1 - ^ 2 ))(1 - ^ 2 ) ( ^ 1)/2 5i ^(l - * 2 )" ) ^ 



2 2i/ r(z/ + i) 2 



= / fl _ z 2 )(^-l)/2 f _S (l _ 2X1/ 

2 2 -r(z/ + i) 2 7 1 j Vd^ -1 

^ 2 x d 2 d (/i-l) 2 \ 

x((l-, 2 )(^)/ 2 ^(l-,Y)dz 

= (^+i)- M ( M -i)) y pr\zfdz 

where we used the differential equation 



(1 - ^)-^ - 2,^^ + + 1) - Yf^J = (B.16) 

obeyed by the Legendre functions. Thus, 

/ P^fdz = 1 / Pj*\z?&z 

J (u-n)(y + n + l) J 

and starting from 

r P- V (z) 2 dz 



T(2v + 2) 

which follows from (jB.llj) the proof by induction is complete. □ 
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